Dependence of the MHD shock thickness on the finite electrical conductivity 
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The results of MHD plane shock waves with infinite electrical conductivity are generalized for 
CO ' a plasma with a finite conductivity. We derive the adiabatic curves that describe the evolution of 

, the shocked gas as well as the change in the entropy density. For a parallel shock (i.e., in which 

' the magnetic field is parallel to the normal to the shock front) we find an expression for the shock 

thickness which is a function of the ambient magnetic field and of the finite electrical conductivity of 
Q_) , the plasma. We give numerical estimates of the physical parameters for which the shock thickness is 

O ' of the order of, or greater than, the mean free path of the plasma particles in a strongly magnetized 

I plasma. 
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I. INTRODUCTION 



Previously, shock waves in plasmas, both relativistic and non-relativistic, were studied assuming ideal MHD, [1-3]. 
Although this theory is suitable for studying most astrophysical shock waves, such as those in hot rarified astrophysical 
plasmas, where the electrical conductivity is extremely high and the magnetic field is weak, it is interesting to study 
^ ' the effect of the simplest dissipative process in non-ideal MHD, that due to a finite value of the electrical conductivity 
^ ' a. We generalize here the results for planar shock waves to non-ideal, non-relativistic MHD. The junction conditions 
Q_i that must be satisfied across a shock wave in a non-ideal plasma are given in Sec. H. We then derive the adiabatic 
^ [ curve with corrections due to a and find the corresponding change in the entropy density across the shock. We 
O ' perform this analysis for an oblique shock (i.e., in which B is neither parallel nor perpendicular to the normal to the 
c/3 , shock surface) as well as for a parallel one (i.e., in which B is parallel to the normal to the shock surface). In Sec. HI, 
we find a closed expression for the shock thickness for a parallel shock in a strongly magnetized plasma, and estimate 
its value for some physical situations. Our conclusions are discussed in Sec. IV. 



II. JUNCTION CONDITIONS AND ADIABATIC CURVE 

For non ideal MHD, the fluxes of mass, energy, and momentum are given by Eqs. (1-4), respectively, [1,4], 
(N ■ pv, (1) 
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Tl,k^ pViVk+pS^k- ^\B,Bk- -B^Sik] , (3) 



q = ptJ ( -w^ + w ) + —B X [vx B] - ^^_y_ B x ( V x S ) , (2) 

and 



where p is the fluid density, v the velocity, w the enthalpy per unit mass, B the ambient magnetic field, and p is the 
fluid pressure. The electric field is 
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{v X B^ X B, (4) 
where Ohm's law in its simplest form [2,5] was used. 



A. Junction conditions 

We assume a two dimensional, planar, shock wave in the y — z plane. The normal to the transition surface is in the 
-X direction. The velocity field can be decomposed into perpendicular and tangential components to the surface of 
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transition, v = (u^, vt). It is assumed that all quantities vary as a function of x . Let be a unit vector normal to 
the transition surface. We then have the hydrodynamical junction conditions, [4], 



[pt;^n^] = 



(5) 



[q.n^] = 



(6) 



(7) 



where i is a; (t) [normal (tangential) to the shock surface] and [] means the difference between the value of the 
corresponding quantities far upstream (which wc denote by subscript "1") and the value at some point in the shock 
(no subscript). It is also assumed that both far upstream and far downstream, all gradients vanish (i.e., the fields and 
flows are uniform). 

In the MHD case that we arc considering, Eqs. (5), (6), and (7) must be suplemcntcd with the electromagnetic 
junction conditions, i.e., that the normal component of the magnetic field and the tangential component of the electric 
field must be constant across the shock surface: 



[Bn] = 0, 



(8) 
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Eq. (5) states that the mass flux along x is conserved, i.e., pv^ = j = const. We write p = 1/V, where V is the 
specific volume, and replace Vx = jV in the other junction conditions, obtaining 



^ ■2ir2 , ^ 2 . 
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Bt 



= 0, 
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B. Adiabatic curve and entropy density change 

We derive the expression for the adiabatic curve and the corresponding entropy density change. Prom Eqs. (12) 
and (13), we obtain 
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Prom Eq. (12), we have [vt] = Bn Bt /A-kj. We can therefore complete the squares in Eq. (10 ), obtaining 
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Prom Eq. (12), the third term of Eq. (15) is zero and we are left with 
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Using Eq. (14), we can write the third term of Eq. (16) as 

' $[Bn = l[ym+iiy-y^)B-B^ 
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Prom momentum conservation, we obtain 
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Using Eqs. (17) and (18) in Eq. (16) and w = £ + pV, we have 

e-ei + l{p + Pi){V- Vi) + ^(^Bt- Buf {V - V^) 



- (V X 5) . (St + Bu) . 



(16) 
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The first four terms are found in the equation for ideal MHD (e.g. Ref. [1]), while the last two are due to the 
finite electrical conductivity. The first term in the square brackets is due to the fact that B„ ^ 0; the second is the 
contribution from the tangential component of the magnetic field. 

To obtain the entropy density change, we follow the procedure found in the standard literature, [4] and develop 
V — Vi in powers of {p — pi). We also expand (w — wi) in powers of {p — pi) and to first order in powers of (s — Si). 
The resulting expression is 
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The first two terms are found in the equation for ideal MHD and the last two are the corrections due to a finite <j. 

In the following section, wo repeat the calculations for a perpendicular shock and find that in the expressions for the 
adiabatic curve and the entropy density change, the only term present which depends on the conductivity is the last 
one. 

C. Perpendicular shock 

Shock waves in a plasma permeated with a magnetic field show several features, of which the most well known is 
related to the orientation of the magnetic field with respect to the shock plane. Although perpendicular shocks can be 

considered to be a special case of oblique shocks, it is interesting to write the simplified expressions for the junction 
conditions explicitly, and re-derive the adiabatic curve and the entropy density change for this case. 



1 . Hydrodynamical and electromagnetic junction conditions 

Por perpendicular shocks Bn = 0, so that the junction conditions now read 

[pva;] = 0, 
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= 0, 



[pvtVx] = 0^ vt=vn, 
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2. Adiabatic curve and entropy density change 
Proceeding as for an oblique shock and defining 
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we obtain 
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There is now only one term that depends on the dissipative properties of the plasma, while for the oblique case, we 
had two such terms. The missing term is related to the normal component of the magnetic field. 



III. THICKNESS OF THE SHOCK WAVE 



The calculation of a general expression for the shock thickness is very difficult, if not impossible. However for a 
perpendicular shock it is possible to calculate the shock thickness exactly. We then have Bn = 0, vt = and consider 
a coordinate system in which the only non-zero component of the magnetic field is By = B, [1]. In this case, the 
equation V.B = is satisfied identically. The unidimensional ideal MHD equations are 



(29) 



dp d , . „ 

^ + ^ M = 0, 
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From the first two equations, it is easy to see that the ratio B / p = (3 satisfies the equation dj3/dt + Vxdj3/dx = or 
dp/dt ~ [1]. Hence, if the fluid is homogeneous at some initial instant, so that /? = const., then it will remain so at 
all subsequent times. Substituting B = pp in the third equation, we obtain 
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dt 
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Thus, the magnetic field has been eliminated from the equations. The equation for the velocity field, Eq. (32), is 
formally identical to that for the ideal fluid case, provided we define the 'true pressure' as p* = p + P^p^/Sn. We can 
now proceed to evaluate the thickness, following Ref. [4]. We write 
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where -y*^ = {dp* /dp)^ and L* is the damping length, given by Eq. (61) of the Appendix, 
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From Ref. [4] we have 
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Equation (33) can be solved using the procedure in Ref. [4] , obtaining the thickness of the shock wave as 

4cL* 



6* = 



(33) 



(34) 
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where P2 and pi arc the 'true pressures' far downstream and far upstream respectively. Equation (36) is quantitatively 
valid for sufficiently small differences (p2 — Pi)- However we can use it qualitatively to estimate the order of magnitude 
of the thickness in cases where the difference {p2 — Pi) is of the same order of magnitude as P2 and p* themselves. 
The velocity of sound in the gas, Vg (not v*) is of the same order of magnitude as the thermal velocity v. Let A be 
the mean free path of the atoms in the plasma. Then from dimensional analysis, the electric conductivity can be 
estimated as, a ~ jv/X ~ jVs/X, where 7 takes into account anomalous effects and can have a value 10~^ < 7 < 1 ^• 

'03/ 

relations in Eq. (36) we obtain 



In Eq. (34), we take B^^/ {B^^ + 4it pov"^) ~ B^y/p* and in Eq. (35) «; (Ps -Pi) ~ vf ~ p*/p. Using these 
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The shock thickness is larger than the mean free path when c p B^y > ^p* p. Let us first assume that Bq^ > p. 



We then have p ~ B^y and p < B^y < &(? j^p. As a specific numerical example, consider p 10 gr/cm^ and 
T ^ 10^ K (characteristic parameters at the center of a massive star before collapse). Assuming hydrogen gas. we 
have n ~ 10^^ cm""^ and p ~ nT ~ lO^^erg/cm"^ (for iron nuclei, the pressure would be two orders of magnitude 
smaller). For the above densities and pressures, the shock thickness is larger than the mean free path of the particles 
if the magnetic field is in the interval 10^ G < Boy < 7-^/^10!'' G. If we now assume that B'^y < p, the shock thickness 
is larger than the mean free path if jp^/c^p'^ < B^y. For the above parameters we have 7^/^10^ G < B^y < 10^° G 
^. If neither of the two conditions above are fulfilled, the shock thickness is smaller than the mean free path. This 
means that the MHD approach breaks down and kinetic theory is needed to study the structure of the shock. 



IV. CONCLUSIONS 



In this article, we extended the results of shock waves treated in ideal MHD to the non-ideal case, in which the 

electrical conductivity is finite. Wc considered Ohm's law in its simplest form, [5], but took into accoimt phcnomeno- 
logically (through the parameter 7) plasma effects that can modify the classical Spitzer electrical conductivity (e.g., 



^We know that in accretion disks, protostars, galactic nuclei and neutron X-ray sources, for example, the plasma cannot 
have ideal Spitzer values for the conductivity and viscosity in order to obtain the observed accretion rates. Therefore it is 
generally assumed that these quantities are highly anomalous (due to turbulence, for example). Another example where the 
assumption of anomalous resistivity is used is in the treatment of solar flares, which arc generally assumed to bo duo to magnetic 
reconnection. If ideal Spitzer values are used for the plasma in solar flares, reconnection times are ~ lO'' times longer than the 
observed time scales for the flares. In general, plasmas near shocks are expected to be highly anomalous (i.e., 7 <C 1) due to 
turbulence. 

^For 7 ^ 10~®, the magnetic field range is 10^ G < Boy < 10^" G. A magnetic field B > 10^ G can be easily be present at the 
center of a massive star. The magnetic field increases in the collapse of the core of a massive star to B > 10® (for the collapse 
to a white dwarf) and to B > 10^ (for the collapse to a neutron star). Thus, at the center of a massive star, we may expect 
that the magnetic field varie from 10^ G to 10^° G during the collapse of its core and the start of a supernova explosion. 
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turbulence). The expressions for the adiabatic curve and the entropy density change across the shock were generalized. 
Finally, we derived the expression for the shock thickness for a finite conductivity in the case of a parallel shock in 
strongly magnetized plasmas. The conditions that the ambient magnetic field must satisfy for the thickness to be 
of the order of the particle mean free path were estimated. We found that these conditions can be fulfilled for the 
plasma expected in the origin of a supernova explosion, [6,7]. Extensions of the results presented in this paper, using 
a more general Ohm's law, as well as to relativistic shocks, are presently under investigation. 
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V. APPENDIX 

In this appendix, we sketch the derivation of the expression for L* , the damping length used to calculate the shock 
thickness. Neglecting the displacement current (which is a good approximation in non-relativistic electrodynamics), 
the evolution equation for the magnetic field in a medium with electrical conductivity a moving with a velocity v is 

Adding the equations for the fluid, which we assume has neither viscosity nor thermal conduction, we have 

1^ + V. ipv) = 0, (39) 



^+(v.v)v = --yp--^Bx(vxB). (40) 
at \ J p A-Kp \ J 



A. Hydromagnetic waves 

Let us assume that B = Bo + b, p = pq-\- 5p, p = Pq-{- 5p, and v = 5v. Replacing these terms in the above equations, 
keeping only terms to first order in the perturbations, expanding the density in powers of the perturbation in the 
pressure (i.e., Sp = Sp/v"^ + [d'^ p/ dp^) ^Sp^ , where Vs is the sound velocity of the medium) and taking the Fourier 
transform of the equations, we obtain 

-uj5po + v^pokMo = 0, (41) 



— (lo + i-^ — P ) 6o = fc X 



47rcr 



5vq X Bq 



(42) 



-LoSvo = Spo - Bo X (kx ba . (43) 

po Anpo V / 

From Eq. (41), we find 5po = v^po (^k.Svo^ /u) and using this in Eq. (43 ), we obtain 

(44) 

We define the scalar phase velocity as m = Lv/k, assuming that k is along the x-axis, (i.e., k = kx) and that Bq is 
in the x — y plane. Writing the previous equations in its components, we have 
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^po = ^5t;ox, (45) 
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v'^\ 1 

u ] Svox = -. boyBoy, (46) 



u ) 47rpo 



u^VQy = - j^boyBox, (47) 
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M + ^'^^^ ) ^oy = ^'^oxBoy - SvoyBox, (48) 



u^wo^ = --^ftoz-Boa;, (49) 



c2 



w + ) boz = -SvozBox- (50) 



B. Generalized Alfven waves 



Using Eqs. (49) and (50) we obtain the compatibility relationship 

. C A/ Bri^ 

I n 

Ana 47rpo 

from which we obtain 



u^+i^u-P^^O, (51) 



1 BL c^k^ 



c 



""-^2^/^ 16^^ 'i^- ^^^^ 

From Eq. (52), the phase velocity is a complex number if a is finite. Rewriting ua in terms of iv, we obtain the 
dispersion relationship: 



W = - \ 1 - ^1:75-— K - I — 



Wo take the plus sign in cq. (52) since the frequency is a positive quantity. The fact that the imaginary part is 
non-linear in k means that the Alfven waves are damped and dissipated as a function of k. For a ^ oo we recover 
the known dispersion relationship for ideal MHD. 

Assuming that the second term in the square root in Eq. (53) is much smaller than unity, the group velocity is 

du! Box .(?k 

When the electrical conductivity is infinite, we recover the known ideal MHD result, vai = Bqx/'^s/i^Po- 
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C. Generalized magnetosonic waves 



Prom Eqs. (46), (47) and (48), we obtain the generalized dispersion relationship for magnetosonic waves, 

Anpo J 47rpo 47rcr Ana 



(55) 



This relationship can be inverted to obtain ui = u{k) . However, for our purposes, it suffices to consider the dissipative 
terms as a correction to the ideal dispersion relationship. 



vi ± 
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k^^vl,k\ 



(56) 



The plus sign corresponds to fast magnetosonic waves, while the minus sign to slow magnetosonic waves. Replacing 
LOi in the last two terms of Eq. (55), we obtain 
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If the electrical conductivity is large (but not infinite), wc have 



2 2 ^ * C WgO 
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where v%q = (Bq/Attpq + v^) and v%q^ = Bq^/ttpo- We thus have 

u ~ Vgok — icL*k^, 

with 

I "bo "BOx"s> 

where L* is the damping length. For a perpendicular shock, (i.e., i?ox = 0), we have 



L* = 



Bit 



8na {Bit + ^t^Pov^) 



(57) 
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